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Abstract: To delay resistance development to Bacillus thuringiensis (Bt) plants expressing 
their own insecticide, the application of the Insect Resistance Management strategy called 
"High Dose/Refuge Strategy" (HD/R) is recommended by the US Environmental 
Protection Agency (US EPA). This strategy was developed for Bt plants expressing one 
toxin. Presently, however, new Bt plants that simultaneously express two toxins are on the 
market. We used a mathematical model to evaluate the efficiency of the HD/R strategy for 
both these Bt toxins. As the current two-toxin Bt plants do not express two new Cry toxins 
but reuse one toxin already in use with a one-toxin plant, we estimated the spread of 
resistance when the resistance alleles are not rare. This study assesses: (i) whether the two 
toxins have to be present in high concentration, and (ii) the impact of the relative size of 
the refuge zone on the evolution of resistance and population density. We concluded that 
for Bt plants expressing one toxin, a high concentration is an essential condition for 
resistance management. For the pyramided Bt plants, one toxin could be expressed at a low 
titer if the two toxins are used for the first time, and a small refuge zone is acceptable. 

Keywords: insect resistance management; population dynamics; population genetics; 
Bacillus thuringiensis; High dose/Refuge strategy; Bt plants; two toxins 
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1. Introduction 

Bacillus thuringiensis (Bt) was discovered in 1901 in Japan in diseased silkworms and has been 
produced on an industrial scale for its insecticidal properties since 1938; however, only in 1953 was it 
shown that its insecticidal properties were the result of the parasporal inclusion of a crystal protein 
called Cry toxin [1]. Unlike many synthetic insecticides which traverse the integument and act on the 
insect's nervous system, these toxins have to be ingested [2]. The proteolytically activated Cry toxin 
binds to specific receptors on the brush border membrane of the midgut. Pores are formed and the 
affected cells swell and are lysed. The death of the insect follows within hours or days [3]. These 
insecticidal proteins are active ingredients of some environmentally friendly insecticides because they 
are harmless to natural enemies and other nontarget organisms due to their narrow host specificity and 
short life in the environment [4]. Today, Bacillus thuringiensis is considered the most successful insect 
pathogen and presently constitutes about 2% of the total insecticidal market [5]. 

A Cry gene was first introduced into tobacco and tomato plants in 1987. These transgenic hybrids 
expressing the insecticidal toxin (Bt plants) are able to withstand insect attack due to exceptional 
expression levels of a Cry toxin [6,7]. Introduced on the market in 1996, Bt plants expressing one toxin 
covered about 50 million hectares (ha) worldwide in 2010 and have been planted on >200 million 
hectares (ha) since 1996. More precisely, Bt corn and Bt cotton covered 42 million ha in 2007 and their 
target insect are Lepidoptera and Coleoptera [8]. One risk linked to this extensive use is an increasing 
selection of resistant insects to Bt plants, which could limit the use of Bt plants in time and space [9]. 

To delay and manage resistance outbreaks in the field environment, a frequent solution is to 
develop and apply Insect Resistance Management (IRM) strategies. These preventive strategies have 
to be applied before resistance spreads widely in the population. Numerous IRM strategies have been 
proposed: (1) mixtures, mosaics or rotations of transgenic plants; (2) time- or tissue-specific toxin 
expression; (3) low doses of toxins in combination with natural enemies; (4) co-expression of different 
Cry genes; (5) high expression of the toxin (dose) with refuge of non-i^ plants [9]. This last strategy, 
the so-called "High Dose/Refuge Strategy" (HD/R), requires planting a "refuge zone" composed of 
non-Bt plants suitable for the target pest and in close proximity to a "Bt zone" producing a high 
concentration of Cry toxin (25 to 50 times the LD99) [10]. 

The HD/R strategy assumes that resistance is an autosomal trait controlled by a single locus with 
two alleles (the rare resistance allele: Freq Ar < 10 — the susceptible allele: As), with no sex linkage 
and no maternal effect, as largely supported experimentally [11-13]. Susceptible adults emerge from 
the refuge zone and resistant homozygotes (ArAr) from the Bt zone. The assumptions of this strategy 
are that mating is panmictic and egglaying is random between the refuge and the Bt zones. If 
inheritance of resistance is functionally recessive, heterozygous progeny produced by such matings 
will be killed by ingestion of high concentration toxin in the Bt zone at the same rate as homozygous 
susceptible larvae. By killing heterozygotes the high dose crop can greatly slow down the increase of 
the resistance allele frequency. 

A substantial number of studies have evaluated the impact of these strategies on resistance 
development. A key tool in this field is mathematical or computer modelling, and various models have 
played an essential role in building a conceptual framework for resistance management [14-19]. 
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Today, Bt cotton and Bt maize still dominate the market and the new hybrids registered in the 
United States collectively produce 18 different combinations of 11 Bt toxins [20]. The second- 
generation Bt plants simultaneously produce two dissimilar insecticidal toxins active against the same 
pest [21]. This "pyramidal" approach has been on the market since 2002. US EPA has registered 
Bollgard II cotton (CrylAc in combination with Cry2Ab2; Monsanto) in 2002, Widestrike cotton 
(Cry IF and CrylAc; Dow Agro-sciences) in 2004 and, MON 810 x TCI 507 Corn (CrylAb and 
Cry IF; Pioneer/Dupont) in 2010 [21]. The dual-toxins plants provide a substantially better control of 
insect pests [22]. They should delay the spread of resistance when selection for resistance to one of the 
toxins does not cause crossresistance to the other toxin because resistance to each toxin is initially rare 
and cases of double resistance will be extremely infrequent. This was demonstrated by theoretical 
models and by experiments [23-31]. However, IRM strategies remain essential because dual resistance 
should be triggered in large populations by uniform and long-term simultaneous exposure to the two 
toxins over large areas [32]. The widespread adoption of these transgenic varieties has thus greatly 
increased the opportunity for resistance selection in the field [33]. 

Initially, the HD/R strategy was developed for one-toxin Bt plants. This study evaluates its 
efficiency for pyramided Bt plants. 

In Australia, Bollgard® was directly replaced in 2004/2005 with Bollgard II®. In the USA, during 
the transition to the second-generation Bollgard II®, both the single- and two-toxin Bt cottons have 
been grown since 2003. The process ended in 2010, since the US EPA registration for Bollgard I® 
cotton expired at the end of 2009. However in India, Bollgard® and Bollgard II® coexists since 
2006 [34,35]. Zhao and colleagues [30] studied the impact of the concurrent use of Bt plants 
expressing either one single, or two Cry toxins, and sharing one of these. They concluded that this 
procedure selects for resistance to two-toxin Bt plants more rapidly than the use of two one-toxin Bt 
plants alone. Therefore, two toxin Bt plants should not be deployed simultaneously with one-toxin Bt 
plants if they share one similar Cry toxin. With our model, we evaluate the spread of resistance when a 
two-toxin Bt crop is planted alone after a one-toxin Bt crop (sharing a similar Cry toxin). In this case, 
the resistance allele frequency to the shared toxin should not be rare (>10~ 3 ) if the resistance is 
developed for the toxin already used [12]. 

One essential condition of the HD/R strategy is the high toxin concentration in the Bt plant. In this 
study, we assess if this condition is necessary for the pyramided Bt plants.We also measure the impact 
of the relative size of the refuge zone on changes in resistance and in the population density in our 
study area. 

For these purposes, a deterministic and discrete time model in two parts was developed. The first 
part is based on population genetics theory and allows us to follow the frequency of resistance at each 
generation. The second is a population dynamics model, which simulates population changes. This last 
part is essential to assess the impact of an IRM strategy on pest population density because controlling 
the pest populations does, in itself, limit crop damages. 
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2. Materials and Methods 

2.1. Parameters 

The effectiveness of any IRM strategy is related to two kinds of parameters: the operational 
parameters characterizing the strategy and the biological parameters characterizing the target 
insect [36]. The operational parameters include the relative size of the refuge zone (v) in relation to the 
whole field and the selection pressure exerted by the Cry toxin on susceptible insects (slfrA, sBtB). 
The biological parameters are the initial frequency of the resistance allele for each locus (ArFreq, 
BrFreq) and its dominance (hAr, hBr), the fitness cost associated with resistance (fcostA, fcostB) 
and its dominance (hfcA, hfcB), the initial population size (nzero), the population intrinsic growth rate 
(r) and the field carrying capacity (K). Each parameter is associated with a default value taken from 
the literature and measured on a lepidopteran insect pest (see Table 1). 

Table 1. Operational and biological parameters introduced in the simulation model: their 
symbol, default and tested values. 



Symbol Default-Value Tested-Value 







Refuge zone proportion 


V 


0.05 [37] 


0.05-0.10-0.20- 


-0.30-0.40 


a 


U 

o 


Selection of the toxin A 


sBtA 


1 [37] 


1-0.93-0.50 




es 
■- 


'aramc 


Selection of the toxin B 


sBtB 


1 [37] 


1-0.93-0.50 




a> 

a 


Field area (hectare) 




260 [38] 


na 




O 




Plants/hectare 




67,000 [38] 


na 








Initial Ar frequency 


ArFreq 


1.5 x 10~ 3 [39] 


1.5 x 10~ 2 -1.5 x 


io- 1 






Initial Br frequency 


BrFreq 


1.5 x 10 3 [39] 


1.5 x 10 2 -1.5 x 


io- 1 


■- 




Ar dominance 


hAr 


0 [12,40] 


0-0.23-0.53 




9) 
= 




Br dominance 


hBr 


0 [12,40] 


0-0.23-0.53 




1 Parai 




Fitness cost associated to Ar 


fcostA 


0.15 [41] 


na 






Fitness cost associated to Br 


fcostB 


0.15 [41] 


na 




)logica 




Fitness cost dominance associated to Ar 


hfcA 


0 [42-44] 


na 






Fitness cost dominance associated to Br 


hfcB 


0 [42-44] 


na 




s 




Initial individual number/ha 


nzero 


50,000 [38] 


na 








Intrinsic growth rate 


r 


0.15 [45,46] 


na 








Carrying capacity/plant 


K 


22 [38] 


na 





2.1.1. Resistance Alleles 

There is one single locus associated with resistance for each of the two toxins (toxin A and toxin 
B). Each locus is characterized by two alleles: the resistant and the susceptible allele (Ar-As and 
Br-Bs) and we assumed that there was no linkage. 

Among the large number of insecticidal Cry many exhibit cross-resistance (so, resistance to one 
toxin confers an increased resistance to a second Cry toxin [47]). Heliothis virescens (Lepidoptera: 
Noctuidae) (Fabricius) individuals selected with one single type of Cry toxin also displayed a 
broad-spectrum Bt resistance. There was a non-surprising crossresistance between CrylAc and 
CrylAb (these toxins are similar in structure and toxicity) and an unexpected crossresistance between 
CrylAc and Cry2A (the amino acid sequences are very different and the toxins differ in their mode of 
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action) [48]. Crossresistance between CrylAb and CrylAc is also recorded for Plodia interpunctella 
(Lepidoptera: Pyralidae) (Hiibner) [49]. Generally, complete crossresistance (selection with one toxin 
leads to an identical change in resistance levels of a second toxin) is possible but in many instances, 
selection with a toxin only confers partial crossresistance to other toxins [48]. The combination of the 
toxins has to be approached very carefully for an appropriate resistance management [28]. In this 
study, we do not consider the crossresistance 

As an initial condition, the resistance allele frequency is low (ArFreq = BrFreq: 0.0015) but the 
efficiency of the HD/R strategy is evaluated when resistance is not rare (ArFreq, BrFreq: 0.015 
or 0.15). 

2.1.2. Relative Size of the Refuge Zone (v) 

Theoretically, a refuge zone should produce 500 susceptible insects for every resistant insect that 
survives on a Bt crop [50]. Determining the optimally efficient size of the refuge zones, their 
arrangement and their acceptability for the users or sellers of Bt crops is, today, a matter of some 
debate [51]. CrylAc cotton growers of the USA have to choose among different structural options: 5% 
external and unsprayed refuge; 20% external and sprayed refuge and 5% embedded refuge [52]. On the 
other hand, Monsanto Company has petitioned the US EPA to eliminate the refuge where Bollgardll® 
cotton is planted, as natural refuge is assumed enough to maintain susceptible insects [53]. The US 
EPA has approved this alternative [54]. In our study, we only consider the externally-structured and 
unsprayed refuge v with a default value of 5% in our standard conditions but other values (10%, 20%, 
30% and 40%) are tested. 

2.1.3. Mortality Exerted by the Bt Plants on Susceptible Insects and Resistance Allele Dominance 

According to the initial assumptions of HD/R strategy, Bt plants produce high toxin titer (sBtA = 1) 
and the resistance allele is recessive (hAr = 0). This will eliminate the heterozygous offspring and slow 
down the increase of the Ar frequency [9], because when an allele is rare, it occurs almost exclusively 
in the heterozygous state. Therefore, during early generations, the resistance allele dominance is the 
main determinant of the rate of evolution [36]. Basically, the HD/R strategy requires a resistance 
controlled by a recessive single locus with two alleles (As-Ar) [11-13]. In most cases studied so far, 
the inheritance of high levels of resistance to Bt toxins in field- and laboratory- selected populations is 
recessive [55,56]. Although the general consensus is that only the operational factors can be 
manipulated [36], experimental data show relationships between Bt toxin concentration, larval 
mortality and functional dominance [11,40,57,58]. Dominance of resistance is not a fixed value and 
presents a plastic response depending on environmental parameters: the recessivity of resistance is 
associated with more demanding environments [59]. Resistance to CrylAc in Pectinophora 
gossypiella (Lepidoptera: Gelechiidae) (Saunders) is codominant at low concentrations (sBtA = 0.5; 
hAr = 0.53), partially recessive at middle concentrations (sBtA = 0.93; hAr = 0.23) and completely 
recessive at high concentrations (sBtA = 1, hAr = 0) [13]. These three situations are evaluated in this 
study to assess the impact of the toxin concentration when Bt plants express two toxins 
simultaneously: HsBtA or sBtB: 1; hAr or hBr: 0 (recessive), 2/sBtA or sBtB: 0.93; hAr or hBr: 0.23 
(partially recessive) and 3/sBtA or sBtB: 0.50; hAr or hBr: 0.53 (codominant). 
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2.1.4. Fitness Cost of Resistance and Fitness Cost Dominance 

Selection for any novel trait can induce changes correlated with other characters, mainly through 
the ability of a single gene to affect more than one trait (pleiotropy) [60]. This antagonistic trade-off in 
fitness caused by the presence of a resistance gene will be referred as a "fitness cost." Here, we 
suppose that it only affects survival. The fitness cost and its dominance are central for predicting the 
behavior of the resistance gene because the cost increases the advantage of SS over RS individuals in 
refuges and favors the decline of resistance [61]. The default values of these two parameters are the 
same as in simulations of Tabashnik et al. [42] and of Vacher et al. [62]: fcostA = fcostB = 0.15 and 
hfcA = hfcB = 0. 

2.1.5. Fitness Estimation 

The fitness is the x ability to survive and reproduce, and thus have offspring in future 
generations [47]. Its value is associated with the genotype (ArAr, ArAs, AsAs for toxin A or BrBr, 
BrBs, BsBs for toxin B) and to the environment. Fitness is calculated for each genotype in the Bt or 
refuge zone following the Lenormand equations (equation 1) [63]. 

w ifcfr = l-/cost4 

w ArAs = 1 - hfcA xfcostA-(\-hAr)x sBtA xgA(x) ( i ) 

w AsAs =l-sBtAxgA(x) 

The function gA(x) equals one if the insect is on Bt plants expressing the A toxin and equals zero 
otherwise. The parameter hAr is the dominance level of the resistance allele associated with toxin 
selection [62]. Similar considerations are made for the B toxin. 

As there is no crossresistance between the two toxins and no pleiotropic effect, the fitness of the 
genotype of the locus A is independent of the fitness of the genotype of the locus B. Thus, with 
two-toxin Bt plants, the global fitness is the product of the fitness of both genotypes (e.g.: 
w ArAr BrBr = w ArAr x w BrBr). 

2.1.6. Initial Individual Number (Nzero); Carrying Capacity (K); Population Intrinsic Growth Rate (r): 

For the values of the parameters characterizing the population dynamics in this system, we use the 
published data concerning the main Lepidopteran pests (see Table 1). K and r are independent of the 
zone in the field {Bt or refuge zone), and of the genotype [45,46]. 

2.2. Model Description 

The considered system is a section of a homogeneous region of 260 ha, as proposed by 
Guse et al. [38]. This field is closed and composed of two adjacent zones: the Bt zone and the refuge 
zone. The insect population is uniformly distributed among the zones and cannot leave the area. Each 
generation is divided into a succession of simple stages and the generations are discrete and continuous 
(i.e., no diapause) (see Figure 1). The time unit is the non-overlapping generation. 
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Figure 1. Stages composing the model. 



0. Panmictic mating at the field scale 



1. Random egglaying at the field scale 














2. Larval selection by the refuge zone 






2. Larval selection by the Bt zone 



3. Situation at the end of the generation 



At step 0, there is panmictic mating at the field scale followed by random egglaying between the 
two habitats (Step 1). Oviposition choice tests indicate that females are unable to discriminate between 
the two plants varieties and show no oviposition preferences [64,65]. There is no hatching difference 
between the two zones [64]. After eclosion, the larvae start feeding and are subjected to selection 
pressure according to their genotypes and to the zone in which they develop (Step 2). Toxin exposure 
occurs only for larvae. Next, the model calculates the genotype frequencies and the population sizes 
(Step 3). The adults emerge and the cycle starts again at Step 0. Subject to the condition of panmictic 
mating and random egglaying, the genotype frequency is present in Hardy- Weinberg proportions in 
each zone at the beginning of each generation [47]. This Hardy- Weinberg equilibrium is disturbed 
from the first step of selection but will be restored at Step 0 in the next generation. There is an 
equation for each stage, implemented in R [66]. 

Usually, one prerequisite of the refuge theory is that the resistance level of the pool of susceptible 
insects from a refuge must remain unaffected by an outward migration from the Bt zone 
population [67]. Our model is more realistic: the random egglaying of the resistant females allows the 
introduction of resistant offspring in the refuge; thus, the refuge is not 100% susceptible but can be 
"polluted" by resistant alleles. 

The first part of this model is based on population genetics theory and allows us to follow the Ar 
frequency at each generation for one-toxin Bt plant or the resistance haplotype frequency (ArBr) for a 
two-toxin Bt plant. This model is inspired by the model of Mallet and Porter [68] and adapted from the 
equation of Crow and Kimura [69], which indicates the direction of change in the resistance allele 
frequency at each generation (t). 

The second part of the model is based on population dynamics theory and allows us to follow the 
quantitative population changes. Here, reproduction is density-dependant and the population growth 
follows Ricker's formulation of a classic discrete logistic population model (equation 2) [70]. The 
number of pests at time N( t +\) depends on the population size at time N®, on the intrinsic rate of growth 
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(r) and on the carrying capacity (K). 



r* 1- 



N, 



(2) 



2.3. Model Output 

The model investigates the effect of each parameter on the evolution of both resistance and 
population size. Two indicators are used as standards for comparing simulations. 

The most common indicator of resistance spread is the number of generations required to reach 
50% frequency of the resistance allele Ar for Bt plants expressing one toxin and 50% frequency of the 
resistant haplotype ArBr for the pyramided plants (GF50) [36,68]. This is a convenient measure of 
resistance, which is independent of assumptions regarding population growth [31]. 

The second indicator is the percentage of population decrease between the first two generations and 
after ten generations. Such an indicator is faster to obtain: only ten generations are required. Secondly, 
it can be positive or negative, indicating that the population density is increasing or decreasing, 
respectively, making it more easily interpretable. We consider that the impact of the IRM strategy on 
population density must be assessed, because it is an important indicator of the sustainability of 
the strategy. 

3. Results and Discussion 

3.1. Efficiency of the HD/R Strategy for Bt Plants Expressing High Dose of One or Two Toxins, with 
Initially Rare Resistance Alleles 

Our simulation results are consistent with the HD/R strategy assumptions, with published 
theoretical models and experimental data [9,23-27,31,71]. The spread of resistance in the population is 
efficiently delayed by Bt plants producing one toxin with a high concentration level. The critical 
threshold of 50% resistance allele frequency in the population (GF50) is reached at the 48th generation 
(see Figure 2A and Table 2A). According to our results, this strategy is more effective with Bt plants 
expressing simultaneously two toxins because resistance is reversed in the population (see 
Table 2B). Only double homozygotes (ArArBrBr) are resistant and resistance to two toxins is a rare 
event (initially ArArBrBr frequency: 5.06 x 10 and ArAr frequency: 2.25 x 10 ). The resistance 
alleles are mainly present in the heterozygote insects which are eliminated by the toxins. This induces 
a decrease of the resistance allele frequency in time. The control of the pest population is similar for Bt 
plants with one or two toxins: there is a reduction of 94% of the population between the first two 
generations and a elimination of the population after ten generations (100% population decrease 
between the tenth and the first generation (see Table 2A,B). 
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Figure 2. Efficiency of the HD/R strategy with Bt plants expressing one or two toxins 
(initial conditions). These simulations are performed for Bt plants expressing one or two 
toxins with the initial conditions: high toxin concentration (sBtA = sBtB =1), recessive 
resistance allele (hAr = hBr = 0), initially rare resistance alleles (Ar Freq = Br Freq = 
0.0015) and a 5% refuge zone. 
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■ 
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3.2. Efficiency of the HD/R Strategy if Resistance Is Initially Not Rare in the Population 

There is no guarantee that the frequency of the resistance allele is rare. In Arizona, a field 
population of Pectinophora gossypiella (Lepidoptera: Gelechiidae) (Saunders) presented in 1997 an 
estimated mean frequency of a major resistance allele to CrylAc Bt cotton ca. 100 times higher than 
the theoretical threshold [12]. In the southeastern United States, 14 field populations of 
Helicoverpa zea (Lepidoptera: Noctuidae) (Boddie) had more than a hundred-fold higher than 
theoretical resistance to CrylAc Bt cotton, and showed a survival increase on leaves of Bt cotton 
plants [20,72-74]. Similarly, high resistance reports have been made for Spodoptera frugiperda 
(Lepidoptera: Noctuidae) (Smith) to Cry IF Bt corn in Puerto Rico and Busseola fusca (Lepidoptera: 
Noctuidae) (Fuller) to Cryl Ab Bt corn in South Africa [75,76]. 
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Resistance to one-toxin Bt plants can be problematic for the pyramided Bt plant because two-toxin 

® 

Bt plants can have the same Cry toxin as the previous crop of one-toxin Bt plants (Bollgardll cotton 

® 

produced Cry 1 Ac and Cry2Ab derives from BollgardI cotton protected by Cry 1 Ac). We consider this 
last situation in our model, with the complementary hypothesis that, due to previous exposure to one of 
the toxins, the resistance allele to this toxin is not rare. 

For the one-toxin Bt plants, the initial frequency of the resistance allele has a high impact on 
resistance spread in the population (see Table 2A, Figure 3A). When the resistance allele is not rare, 
the efficiency of the resistance control drastically decreases. 48 generations are necessary to reach 
GF50 with an initially rare resistance allele (Ar Freq: 0.0015). This threshold is reached after eight 
generations for a 0.015 Ar frequency and after two generations for a 0.15 Ar frequency. Our results 
validate that the HD/R strategy is a preventive strategy, which is only suitable with initially rare 
resistance alleles [31,32]. The control of pest density is also influenced by the initial allele frequency 
(see Table 2A). After two generations, there is a 94% population decrease with 0.0015 and 0.015 
initial Ar frequencies but a 92% population decrease with a 0.15 initial Ar frequency. However, 
whatever the initial resistance allele frequency, the insect population is eradicated after 
ten generations. 



Table 2. Evolution of resistance and population density. 



A. 






GF50 


% Pop. Decrease 




m pianis »ynmesmng yjnt i oxin 


Gen 2/Genl 


Gen 10/Gen 1 


1. 


Initial Conditions 

(sBtA = 1, hAr = 0, Ar Freq = 0.001, v = 5%) 


48 


94 


1 


2. 


Impact of the Ar 
Frequency 


Ar Freq 0.0015 
Ar Freq 0.015 
ArFreq 0.15 


48 

8 
2 


94 
94 

92 


1 

99.9 
99.1 


3. 


Impact of the Cry 
Concentration 


hAr 0-sBtA 1 

hAr 0.23-sBtA 0.93 

hAr 0.53-sBtA 0.50 


48 
7 

20 


94 

86 
39 


1 

99.9 
98.7 



4. v0.05 

vO.l 

Impact of the Refuge 

v 0.2 

Zone Proportion 

v 0.3 
v0.4 



48 94 1 

96 88 99.9 

213 77 99.9 

373 65 98.9 

<500 54 99.8 
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Table 2. Cont. 



B. 










% pop. decrease 






Bt plants synthesizing two toxins 


GF50 


Gen 2/ 


Gen 10/ 












Genl 


Gen 1 


1. 


Initial Conditions 






>500 


94 


1 


2. 


Impact of the Ar 


ArFreq 0.0015 




>500 


94 


1 




Frequency 


ArFreq 0.015 




>500 


94 


1 




(Br: 0.0015) 


ArFreq 0.15 




>500 


94 


1 


3. 


Impact of the Ar 
and Br Frequencies 


ArFreq0.15-BrFreq0.15 


8 


94 


99.9 


4. 


Impact of the Cry 




hBr 0-sBtB 1 


>500 


94 


1 




concentration 


hAr 0-sBtA 1 


hBr 0.23-sBtB 0.93 


>500 


94 


1 








hBr 0.53-sBtB 0.50 


107 


94 


1 






hAr0.53-sBtA0.50 


hBr 0.23-sBtB 0.93 
hBr 0.53-sBtB 0.50 


22 
25 


90 
67 


99.9 
99.9 






hAr 0.23-sBtA 0.93 


hBr 0.23-sBtB 0.93 


21 


94 


1 


5. 


Impact of the Cry 


hAr 0-sBtA 1 


Ar Freq 0.0015-Br Freq 0.15 


>500 


94 


1 




Concentration and 


hBr0.23-sBtB0.93 


Ar Freq0.15-Br Freq 0.0015 


9 


94 


q 




of the Resistance 


hAr 0-sBtA 1 
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Figure 3. Efficiency of the HD/R strategy with Bt plants expressing one or two toxins 
when the resistance is not rare in the pest population. These simulations are performed for 
Bt plants expressing one or two toxins with a high toxin concentration (sBtA = sBtB = 1), a 
recessive resistance allele (hAr = hBr = 0) and a 5% refuge zone. However, the resistance 
alleles are not initially rare (Ar Freq and Br Freq = 0.0015, 0.015 or 0.15). 
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For the Bt plants expressing two toxins, resistance management remains effective when only one 
resistance allele is not rare in the population (see Figure 3B and Table 2B). The high toxin 
concentration {sBtA = sBtB: 1) eliminates all heterozygotes, whereas the homozygous resistant insects 

2 2 

(ArArBrBr) are very rare even if one resistance allele is not infrequent (0.0015 * 0.15 ). In these 
cases, resistance frequency decreases in the population (so, the GF50 threshold is never reached) and 
there is a 94% population decrease after two generations and an elimination of the pest population after 
ten generations. When the two resistance alleles are not rare (Ar = Br frequency: 0.15), resistance 
spread is rapid (GF50: 8) but the population density is still well controlled (94% and 99.9% population 
decrease after two and ten generations, respectively) (see Table 2B and Figure 3C). Therefore, 
although the initial resistance allele frequency for one- or two-toxin Bt plants does not seem to 
influence population density, it can nonetheless have a catastrophic impact on resistance control if the 
two toxins are not initially rare. 

3.3. Efficiency of the HD/R Strategy IfBt Plants Produce a Low Toxin Concentration (for One or 
Two Toxins) 

The efficiency of the high toxin concentration for one -toxin Bt plants on pest management was 
experimentally demonstrated in Arizona with a large-scale, long-term (ten years) study. The extensive 
use of Bt cotton had caused regional suppression of Pectinophora gossypiella, an ecological cotton 
specialist [77]. Thus, regions with substantial Bt fields imposing high mortality could have negative 
effects on population growth of a specific pest with a narrow host range before resistance occurs 
because their diet would be more affected [9,18]. According to Carriere et al. [61], the high dose limits 
the spectrum of resistance mutation conferring resistance in a particular species. Only few mutations 
could generate sufficient effects to overcome the toxin. The most widely observed mechanism in the 
laboratory and in the field is related to a reduced binding affinity of Cry toxins to midgut membrane 
receptors [78]. 

However, experimental studies have suggested that high Cry toxin concentration significantly 
increases the development duration of several insect pests [13,46,79-83]. The slowing down of insect 
development induces an asynchrony in emergence between the refuge and the Bt zone [29,84-86], 
which disrupts random mating and leads to assortative mating between adults in each of the two 
zones [81,86,87]. The reduction of random mating would disturb the effectiveness of the HD/R 
strategy and some researchers therefore recommend using Bt plants with low toxin concentration. 

We test the efficiency of the refuge strategy associated with Bt plants producing a lower toxin 
concentration (sBtA: 0.93 or 0.50) with a resistance allele that is not completely recessive (hAr: 0.23 
or 0.53, respectively) (see Table 2A and Figure 4A). Our simulations suggest that resistance is no more 
efficiently controlled for Bt plants expressing one toxin than for those expressing two (GF50: 7 and 20, 
respectively). The pest density is not as reduced after two generations (population decrease: 86% and 
39%, respectively) but this difference of population control disappear after ten generations (population 
decrease: 99.9% and 98.7%, respectively). 
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Figure 4. Efficiency of the HD/R strategy with Bt plants expressing one or two toxins with 
low concentration. These simulations are performed for Bt plants expressing one or two 
toxins with initially rare resistance alleles (Ar Freq = Br Freq = 0.0015) and a 5% refuge 
zone. But the toxin concentration varies (sBtA and sBtB = 0.50 or 0.93 or 1) involving a 
variation of the dominance of the resistance allele (hAr and hBr = 0.53 or 0.23 or 0). 
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Concerning the two-toxin Bt plants, we first test the association: one high toxin concentration/one 
low toxin concentration (see Figure 4B,C and Table 2B). For sBtA: IsBtB: 0.93, the GF50 threshold 
is never reached and for sBtA:l-sBtB: 0.50, GF50 is reached after 107 generations. A 94% population 
decrease is recorded for these two situations after two generations. Next, we consider Bt plants 
producing the two toxins at lower concentrations and we conclude that less than 25 generations are 
required to reach GF50 (sBtA = sBtB: 0.53, GF50: 25; sBtA: 0.53-sBtB: 0.93, GF50: 22; 
sBtA = sBtB: 0.93, GF50: 21). After two generations, the reduction in population density is also not so 
effective (sBtA = sBtB: 0.53: 67% population decrease; sBtA: 0.53-sflfB: 0.93: 90% population 
decrease) except for sBtA = sBtB: 0.93 (94% population decrease). However, after ten generations, the 
population control is efficient whatever the initial condition (100% population decrease). Thus, 
resistance is efficiently managed when Bt plants express at least one toxin at high concentration. We 
do not recommend producing Bt plants with two toxins at low concentrations unless there are no 
effective alternative methods (biological control, host plant resistance, etc.). 

3.4. Efficiency of the HD/R Strategy If Resistance Is not Rare in the Population and Two-Toxin Bt 
Plants Produce a Low Toxin Concentration 

As previously suggested, HD/R strategy applied for two-toxin Bt plants is efficient even if the 
second toxin is produced at low concentration (sBtA: 1 and sBtB: 0.93 or 0.50) or if one resistance 
allele is rare in the population (0.015 or 0.15). However, what happens if two-toxin Bt plants produce 
one toxin at low titer when a resistance allele is not rare in the population? 

Resistance is still reversed with sBtA: IsBtB: 0.93 associated to Ar frequency 0.0015 and Br 
frequency 0.15 (see Figure 5A,B and Table 2B). Resistance management strategy is also effective with 
sBtA: IsBtB: 0.53 associated to Ar frequency 0.0015 and Br frequency 0.15 
(GF50: 90). 

Figure 5. Efficiency of the HD/R strategy when the resistance is not rare in the population 
and the two-toxin Bt plants express different levels of toxin concentration. These 
simulations are performed for Bt plants expressing one or two toxins with 5% refuge zone. 
The impact of the toxin concentration (sBtA = sBtB = 1), the variation of the dominance of 
the resistance allele (hAr = hBr = 0) and the initial frequency of the resistance alleles (Ar 
Freq and Br Freq = 0.0015 or 0.015 or 0.15) are tested. 
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However, if the rare resistance allele (0.0015) is not associated with the high toxin concentration 
(sBtA: 1), resistance management is less viable (see Figure 5B and Table 2B). With sBtA: 1-sBtB: 
0.93 associated with Ar frequency 0.15 and Br frequency 0.0015, the threshold GF50 is reached after 
only nine generations. With sBtA: \-sBtB: 0.50 associated with Ar frequency 0.15 and Br frequency 
0.0015, the threshold is reached after 20 generations. In conclusion, for Bt plants expressing one high 
toxin concentration combined with one low toxin concentration, it is risky to reuse a toxin already 
produced by a one -toxin Bt plant. The resistance allele frequency associated with this "old" toxin may 
not be rare. Following our results, the rare resistance allele has to be associated with the high toxin 
concentration for optimal resistance management. Introduction of two new toxins in the plants is 
recommended when low concentrations are used to limit resistance risk. Moreover, Zhao and 
colleagues [30] suggest that regulatory agencies consider deregulating one-toxin Bt plants as soon as 
two-toxin Bt plants sharing one toxin are available. 

Whatever the situation, the population decrease is around 94% after two generations and 100% after 
ten generations, in every case (see Table 2B). 

3.5. Efficiency of the HD/R Strategy in Relation to the Refuge Zone Proportion 

A key factor in the HD/R strategy is the refuge zone. The efficiency of the strategy is improved 
with high refuge zone proportion. The larger is the relative size of the refuge zone, the more insect 
resistance is delayed (see Figure 6 and Table 2A). However, for Bt plants with two toxins at low 
concentration, resistance is efficiently managed with a 30% refuge zone but the pest population is not 
so well controlled after two generations (population decrease: 65% for sBtA = sBtB: 0.93; 63% for 
sBtA: 0.93 and sBtB: 0.53; 45% for sBtA = sBtB: 0.53) (see Figure 6C, Figure 6D and Table 2B). 
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However, after ten generations, we observe an elimination of the pest population whatever the 
refuge size. 

Figure 6. Improvement of the efficie ncy of the HD/R strategy for Bt plants with low toxin 
concentration by increase the refuge zone proportion. These simulations are performed for 
Bt plants expressing one or two toxins with initially rare resistance alleles 
(Ar Freq = Br Freq = 0.0015). The impact of the variation of the toxin concentration 
(sBtA and sBtB = 1 or 0.93 or 0.50), the variation of the dominance of the resistance allele 
(hAr and hBr = 0 or 0.23 or 0.53) and the refuge zone proportion (v = 0.05 or 0.10 or 0.20 
or 0.30 or 0.40) are tested. 
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Figure 6. Cont. 
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Today, US EPA recommendations allow a small refuge size (5%). The critical question is whether 
enough susceptible insects will survive in the refuge to provide an effective source of susceptible 
alleles. As a resistance episode progresses, the refuge zone could be slightly "polluted" by random 
oviposition. The resistance level of the refuge zone may increase and make it inefficient. However, 
according to our results, a 5% refuge seems to be sufficient to control resistance. That could increase 
acceptability of this strategy by growers unwilling to sacrifice a large field area as a refuge. On the 
other hand, it has been shown that immigration of native susceptible insects could also restore the 
susceptibility of a population in the Bt zone [88]). Counting on immigration of susceptible alleles 
instead of using refuge zones would be feasible only if immigration patterns of susceptible insects are 
well known and consistent, and, usually one does not know a priori whether such events would occur. 
Furthermore, extremely high immigration is not compatible with overall pest control objectives (due to 
the damage caused) and if immigrants are not susceptible, they will not delay resistance [42]. 

Even though individuals with double resistance may be very rare initially, a refuge is still necessary 
for two-toxin Bt plants to delay resistance (our results, [31,32]). The refuge zone may be costly in 
terms of yield but the use of two-toxin Bt plants with high dose can reduce the amount of refuge 
required to delay resistance and for extended period in comparison with one -toxin Bt plants [31,89]. 
With a 5% refuge zone, one-toxin Bt plants delay resistance during 48 generations (GF50), whereas 
with two-toxin Bt plants, the delay is longer (GF50: never reached in our model). Therefore, two-toxin 
Bt plants have the potential to reduce the requirements for refuges in successful resistance 
management, which could encourage the industry to develop two-toxin Bt plants [90]. However, 
undersized refuges remain risky, for instance when mortalities of heterozygote insects are lower than 
expected. The more prudent way to deploy transgenic crops remains to keep refuges as large as 
economically feasible [31]. 
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4. Conclusions 

Selection by a regular and frequent use of insecticide inevitably triggers resistance in the target 
insects. The unique reasonable hope to delay or avoid resistance development is to use practices 
reducing the intensity of selection pressure. With the HD/R strategy, the refuges reduce toxin selection 
pressure and the high toxin titer in Bt plants purges each generation of as much of the resistance allele 
as possible. This strategy was developed for Bt plants expressing one toxin. 

Our model indicates that for Bt plants expressing one toxin, a high concentration is an essential 
condition for resistance management. The pyramided Bt plants are more likely to control resistance. 
One of the two toxins could be expressed at a low concentration, particularly if the two toxins are used 
for the first time, in which case resistance to these toxins should not be present in the pest population 
and a small refuge zone would be acceptable. 
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